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Quenched Results on B Mesons with NRQCD 

Presented by A. Ali Khan^* 

^Dept. of Physics and Astronomy, University of Glasgow, Glasgow G12 8QQ 

Preliminary results on spectrum and decay constant of B mesons from a quenched simulation at /? = 6.0 on 
an 16^ X 48 lattice are discussed. The heavy quark has been implemented using both an NRQCD and the static 
formulation, the light one with a clover improved Wilson formulation. Both the NRQCD Hamiltonian and the 
heavy-light current consistently include terms of 0{1/Mq), the inverse heavy quark mass. 
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1. INTRODUCTION 

Calculations of the spectrum of B mesons and 
the decay constant fs are of topical interest in 
lattice gauge theory. Spectrum calculations allow 
to check the viability of the lattice methods by 
comparison to experimental results. 

In B systems the mass of the heavy quark 
is > 1 in lattice units. If a relativistic action 
like Wilson or clover is used naively for both the 
heavy and the light quark one has to simulate at 
smaller quark masses and then extrapolate the 
results to the h mass. The static approxima- 
tion on the other hand, with an infinitely mas- 
sive heavy quark, is expected to yield too a high 
value for fs compared to a 5 meson with re- 
alistic quarks and cannot reproduce heavy mass 
dependent splittings correctly. In NRQCD the 
heavy quark is treated in the first approximation 
as nonrelativistic, and corrections due to relativ- 
ity are introduced systematically order by order. 
This methods allows us in principle to study B 
systems directly, without the need of any extrap- 
olation. For heavy-light systems the parameter 
with respect to which the action is expanded is 
given by the inverse heavy quark mass l/Mq. 

2. SIMULATION 

The nonrelativistic Hamiltonian H used here to 
describe the h quark and the pseudoscalar heavy- 
light current are corrected through order 1/Mq, 
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where Mq is the bare heavy quark mass, at tree 
level. So H looks as follows: 
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where Q is the two component heavy quark 
spinor. The Hamiltonian is tadpole improved, 
i. e. the gauge fields are divided by the fourth 



root of the average plaquette: 



(2) 



Also the heavy-light current is corrected through 
0(1/Mq). This is done by considering the Foldy- 
Wouthuysen transform which relates the small 
components of the heavy quark 4-spinor to the 
non-relativistic spinor Q to the desired order in 
Mq. 
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We extract fs from a matrix element of the tem- 
poral component of the axial vector current be- 
tween a pseudoscalar state i.e. the B meson, and 
the vacuum. This matrix element in NRQCD is 
related as follows to its counterpart in relativistic 
QCD: 



{qi5i4qh\PS)QCD 
1 



{q^Q\PS)NRQCD + 
-{q^-ia-D)Q\PS)NRQCD(4:) 

Q 
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Here, a bare heavy quark mass of Mq = 1.71 in 
lattice units is used. 

This computation was performed on a set of 36 
quenched configuration at /3 = 6.0 on a 16^ x 48 
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lattice, fixed to Coulomb gauge. For light quarks 
we use clover improved Wilson fermions with the 
clover coefficient c = 1, i.e. lattice spacing errors 
removed through 0(a) at tree level, at k values 
of 0.1432 and 0.1440. The critical and strange k 
values for these propagators are Kg = 0.14551(3) 
and Kg = 0.1437(1), the scale determined from 
Mp is a-i = 2.05(6) GeV [1]. Here we would 
like to point out that for quenched configurations 
the scale determination does vary heavily depend- 
ing on the physical quantities considered, as ex- 
pected. The scale obtained from T spectroscopy 
at 13 = 6.0 is 2.4(1) GeV [3]. 

We studied s wave correlation functions in 
pseudoscalar and axial vector channels, both lo- 
cal and smeared at source and sink. The smear- 
ing functions used in our run and at SCRI [2] 
are hydrogen-like wavefunctions (ground and 1st 
excited state for the NRQCD propagators and 
ground state for the static ones). 

3. ANALYSIS AND RESULTS 

3.1. B and Bg masses 

In NRQCD the meson mass m can be calcu- 
lated by adding the nonrelativistic binding enery 
which is extracted from the heavy-light correla- 
tion functions and an energy shift: 



m = A + El 
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The energy shift A contains the renormalized 
quark masses and the zero point of the energy. 
Both can be calculated perturbatively [3]. The 
-E'nrqcd's shown in table 1 have been extracted 
from fits of the ground and excited state smeared- 
local correlators to 2 exponentials. An extrapo- 
lation of the results to Kg and Kg yields 0.437(5) 
and 0.491(6) respectively. Converting this into 
physical units we get a B — Bg mass splitting of 

Table 1 

tlesults in lattice units 



K 


0.1432 


0.1440 


-E'nrqcd 


0.506(6) 


0.482(5) 


AE (B* - B) 


0.021(1) 


0.021(1) 


(static) 


0.238(7) 


0.225(2) 


Zi(NRQCD) 


0.120(5) 


0.113(5) 
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Figure 1. Smeared-local correlation functions 
(crosses: ground state, circles: first excited state) 
with 2-exponential fit from t = 7 to t = 17. 
K = 0.1440. 

0.11(3) GeV, compared to the experimental value 
of 0.096(6) GeV. In Heavy Quark Effective The- 
ory (HQET) this splitting is expected to be inde- 
pendent of Mq . For a further discussion of mass 
splittings in the context of HQET see [2]. For the 
B and Bg masses themselves we get 4.5(4) and 
4.6(4) GeV respectively. This difference with ex- 
periment may be accounted for by the fact that 
we chose Mq in lattice units to give the right T 
mass assuming that the lattice spacing was 2.4 
GeV, the appropriate value for hh systems. The 
dependence of on the scale in the quenched 
approximation means that we have to retune Mq 
for B mesons, where we believe that the light 
scale is more appropriate for most physical phe- 
nomena. 

3.2. B* - B splitting 

To extract the energy splitting AE between the 
^Si(B*) and the '^So(B) state we fit the boot- 
strapped ratio of the smeared-local correlators to 
a single exponential: 



(X Aexpi-AEt). 



(6) 



CsLi'So) 

We obtain AE = 0.021(1) for both k values. 
From HQET one expects this splitting to be pro- 
portional to 1/Mq, so it is not only directly de- 
pendent on but also through the value of 
the heavy mass in lattice units. Having chosen 
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too small a value for Mq for our simulations we 
rescale AE when converting it into physical units, 
which yields a splitting of 37(6) MeV, compared 
to the experimental value 46(1) MeV. If one as- 
sumes that AE is a short distance quantity, one 
might argue that a higher scale is more appro- 
priate here. Using = 2.4(1) GeV one would 
obtain 50(4) MeV in physical units. 

3.3. fB 

The continuum value of fs is related to the 
corresponding lattice quantity through 



(7) 



For sufficiently large times the ratio of the corre- 
lation functions 



Csl/Css ~ Zl/Zs, 



(8) 



So the first step to determine Z^ is to fit the boot- 
strapped ratio of the two correlation functions, 
starting from the time slice where its plateau be- 
gins, to a constant. Zs is obtained from a sin- 
gle exponential fit to the ground state smeared- 
smeared correlator: 



Css ~ Zl exp (-Et) 



(9) 



The results for Z^ in NQRCD and static cases 
are given in table 1. Extrapolating to the critical 
and strange k values gives 0.103(4) and 0.115(5) 
for NRQCD and 0.201(6) and 0.230(7) in the 
static case. In figure 2 we show the ratio of the 
small components contribution to the first term in 
equation 4, which is about 20%. There is no value 
for Zji in NRQCD available yet, whereas for the 
static theory there is a perturbative result of 0.78. 
If we assume that the light scale would be most 
appropriate we obtain in the static approximation 
fs = 0.28(3) GeV and fg, = 0.32(3) GeV, which 
agrees with the UKQCD value [4] within errors, 
as expected. With = 2.4(1) GeV we get the 
substantially higher values of fs = 0.36(3) GeV 
and fB, = 0.41(4) GeV. 

4. CONCLUSIONS 

The good signal for our correlation functions 
and well chosen s state smearing functions enable 
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Figure 2. Ratio R of the small components cor- 
rection to the pseudoscalar correlation function. 
The solid line shows a fit of the ratio to a constant 
from t = 6 to t = 28. 

us to extract NRQCD energies with a good accu- 
racy. The conversion of results to physical units 
suffers from a large uncertainty in the scale; in our 
case the main source of this uncertainty is quench- 
ing. Only given an accurate value for we may 
be able to calculate energies and fs accurately. 
We find that there are large 1 /Mq corrections for 
Zl when going from a static to an NRQCD heavy 
quark where Hamiltonian and axial vector cur- 
rent are corrected through 0(1/ Mq). Once Za 
has been calculated in NRQCD it will be seen 
how this translates into the 1/Mq dependence of 

fB- 
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